August 2008 



Non-compact QED3 at finite temperature: 
the confinement-deconfinement transition 

Roberto Fiore", Pietro Giudice^ and Alessandro Papa" 

° Dipartimento di Fisica, Universita della Calabria 
and INFN, Gruppo coUegato di Cosenza, 
1-87036 Arcavacata di Rende, Cosenza, Italy. 
Dipartimento di Fisica Teorica, Universita di Torino 
and INFN, Sezione di Torino, 
via P. Ciuria 1, 1-10125 Torino, Italy. 



Abstract 

The confinement-deconfinement phase transition is explored by lattice numerical simula- 
tions in non-compact (2-|-l)-dimensional quantum electrodynamics with massive fermions 
at finite temperature. The existence of two phases, one with and the other without con- 
finement of fractional charges, is related to the realization of the Z symmetry. The order 
parameter of this transition can be clearly identified. We show that it is possible to de- 
tect the critical temperature for a given value of the fermion mass, by exploiting suitable 
lattice operators as probes of the Z symmetry. Moreover, the large-distance behavior 
of the correlation of these operators permits to distinguish the phase with Coulomb- 
confinement from the Debye-screened phase. The resulting scenario is compatible with 
the existence of a Berezinsky-Kosterlitz-Thouless transition. Some investigations are pre- 
sented on the possible relation between chiral and deconfinement transitions and on the 
role of "monopoles". 
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1 Introduction 



Non-compact quantum electrodynamics in (2+l)-dimensions (QED3) with fermionic mat- 
ter is a special theory since it plays a role in contexts ranging from condensed matter to 
particle physics. In particular, this theory is relevant for the characterization of the phase 
diagram of high-Tc superconductors in the temperature-doping plane (see, for instance, 
Ref. [1] and references therein). Moreover, it is an interesting theoretical laboratory for 
the investigation of mechanisms of confinement and for the study of the confinement- 
deconfinement transition at finite temperature. 

This paper aims at contributing to clarify some issues related with the latter topic. 
Understanding the mechanism of confinement, identifying the nature of the confinement- 
deconfinement transition and the related order parameter is of central interest in finite- 
temperature non-Abelian quantum field theories, such as Quantum Chromodynamics 
(QCD). Here definite answers are available in the limiting cases of infinite quark masses 
(i.e. in the pure gauge theory) and zero quark masses (i.e. in the chiral limit). Instead, 
many open questions remain for the general case of fermions with finite non-zero mass. In 
non-compact QED3 the symmetry whose breaking determines the transition is the Z sym- 
metry, independently from the fermionic mass, the order parameter being the Polyakov 
loop with fractional charge. Evaluating the vacuum expectation value of the order param- 
eter and its 2-point correlator allows to get, at least in principle, a complete description 
of the phase diagram of the theory for varying fermion masses, coupling constant and 
temperature. 

The only theoretical scenario suggested so far for the phase diagram of non-compact 
QED3 at finite temperature is presented in Refs. [21 E], where it is found that the effective 
action for the temporal component of the gauge field, Ao{x), in the limit of large fermion 
mass becomes equivalent to the sine-Gordon potential. This led the authors of Ref. [2] 
to conclude that there is a phase transition of Berezinsky-Kosterlitz-Thouless (BKT) 
type [4j,[5]. More precisely, there is a critical temperature T^, depending on the fermion 
mass, above which the system is in a deconfined or Debye phase, while below Tc the 
interaction is logarithmic with the distance, i.e. it is Coulomb-like. In this paper we 
intend to verify the above scenario for non-compact QED3 at finite temperature, by 
Monte Carlo numerical simulations, using as probes the Polyakov loop with fractional 
charge and its 2-point correlator. 

Another topic we consider in this paper is the behavior of the chiral condensate across 
the deconfinement transition. In particular, if the chiral condensate exhibits a sharp drop 
when the temperature is increased through the critical value, then we can conclude that 
the dynamical mechanism responsible for deconfinement affects also the chiral symmetry 
of the theory. 

Finally, we present a few numerical results concerning the magnetic monopole density. 
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While monopoles in compact QED3 are undoubtedly related to confinement in the pure 
gauge case [B] and the same is argued also in presence of fermions [3 El E] , in the non- 
compact theory there is no a priori reason for which they could play any role. However, 
following Ref. [lO], a monopole density can be defined on the lattice in the non-compact 
theory exactly as in the compact one and its behavior with temperature can be studied. In 
(3+l)-dimensional non-compact QED, it turns out that there exist a percolation threshold 
for monopole current networks near the chiral transition [TU], thus suggesting a possible 
relation between monopole condensation and chiral symmetry breaking (see Ref. [H] 
for a criticism to this approach). In this paper we present some determinations for the 
monopole density across the deconfinement transition and briefly discuss the possibility 
of a relation between monopoles and confinement in non-compact QED3. 

The paper is organized as follows: in Section [2] we briefly recall the theory in the 
continuum formulation, the origin of the Z symmetry and the conjectured phase diagram; 
in Section[3]we describe the lattice version of the theory and the operators we use as probes 
of the phase transition; in Section H] we present the numerical results and discuss their 
interpretation; in Section [5] we draw our conclusions and sketch the future perspectives. 

2 The continuum QED3 theory 

The continuum Lagrangian density describing QED3 is given in Minkowski metric by 



where F^^, = d^A^ — dyA^ is the field strength, = df^ — igA^ is the covariant derivative, 
g is the coupling constant (or the electric charge) and the fermion fields ipi {i = 1, . . . , Nf) 
are 4-component spinors, defined in such a way that the theory is parity-invarianlll]. The 
relevant information concerning parity and chiral symmetry of this model can be found, 
for instance, in Ref. Section II. 

2.1 The Z symmetry 

In this Section we briefly recall the origin of the Z symmetry which plays a fundamental 
role in this paper. 

The Euclidean partition function of finite temperature QED is invariant under gauge 
transformations [2l, ^3]. For the field Ah{t,x) the gauge transformation 



(1) 



(2) 



together with the periodic boundary conditions in the time direction 



A^{l/T,x) = A,{0,x), 



(3) 



This implies that the fermion is massive, while the photon is massless. 



3 



where T is the temperature, imphes the following condition: 

9,,x(l/T,f)=9^x(0,x). (4) 
Alike, for the fermion field ip{T,x), the gauge transformation 

^'(r,f) =e*5^("'^V(^,^) , (5) 
together with antiperiodic boundary conditions in the time direction 

7/>(l/T,f) = -V'(0,x) , (6) 

implies 

X(l/T,x) =x(0,x) + ^n, (7) 

where n is an integer. Differently, we can say that x{x) is periodic in the time direction 
with period 1/T up to an integer multiple of 2iT/g. In the language of group theory, we 
can say that if G is the group of all gauge transformations and H is the subgroup of 
those gauge transformations which are strictly periodic, then the quotient group G/H is 
isomorphic to Z, the additive group of integers. We refer to this when we say that the 
theory possesses Z symmetry. Note that this is a symmetry of the partition function in 
presence of dynamical electrons. 

In order to study the realization of the Z symmetry, a version of the Polyakov loop 
operator is introduced, whose average is related to the free energy of an external charge 
9- 

ng(f) = e*^^o^^'^^°'^«(^'°'^) . (8) 

Under the action of an element of G/H, the Polyakov loop operator with charge g trans- 
forms as (here we use Eqs.(l2]) and (^^) 

Therefore, in order that the operator defined in Eq. ([H]) be an order parameter, the charge 
g must not be an integer multiple of the basic charge g. One possibility is to follow 
Ref. [12] and study charges that are rational fractions of the fundamental charge, i.e. 
g = g/m, with integer m. As a consequence, Eq. becomes 

= n„(x)e'^ . (10) 



The realization of the Z symmetry tests the ability of the electrodynamic system to 
screen charges which are not integral multiples of the electron charge. Therefore, it is 
related to (fractional) charge screening and confinement in QED in the same way as the 
Zat symmetry does in finite temperature SU{N) pure gauge theory. 

In conclusion, the operator Ilg{x) can be used as an order parameter for confinement in 
Abelian gauge theories even in presence of dynamical charged particles: if the symmetry 
is unbroken, the loop operator averages to zero and the system is in the confining phase, 
otherwise it is in a non-confining phase. 
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2.2 Theoretical expectation for the phase diagram 

In this Subsection, we briefly recall the ideas put forward in Ref. [2J about the phase 
structure of the theory under consideration. 

The authors of Ref. [2] propose the Polyakov loop with fractional charge as order 
parameter and argue that in parity-invariant electrodynamics with fermions of mass M 
there is a deconfinement transition at finite temperature of BKT type. 

The analysis of Ref. p] is based on computing the effective action V{M,gAQ/T) for 
Ao(r, x), which explicitly exhibits the global Z symmetry. By studying this effective action 
they characterize the type of the phase transition. 

At finite temperature, non-compact QED3 contains three parameters (with the dimen- 
sion of mass): the fermion mass M, the gauge coupling (7^ and the temperature T; the 
dimensionless parameter which governs the loop expansion is the smaller between jM 
and c/VT. 

In the large M limit, so that T/M and jM are small with jT finite, they argue that 
the critical behavior of the theory is identical to that of the 2-dimensional sine-Gordon 
potential, which undergoes a BKT phase transition. Therefore they conclude that also 
in QED3 at finite temperature there must be a BKT transition with a critical line in the 
[M/T, (yf^/T] plane, starting at 



(M/T,^7V3^) = (oo,87r), (11) 



from which they obtain 



rpM-^T _ 9_ /-, r,\ 

^crit. -g^- U^J 

They determine also the critical temperature for the BKT transition up to one-loop order: 

T^^[>>T = . (13) 

The second limit considered in Ref. [2] is the high-temperature limit, T ^ M^g'^. The 
analysis of Ref. [2] shows that the Z symmetry is spontaneously broken, which means the 
system is in the deconfined phase. 

In a subsequent paper [3] the same authors characterize better the nature of the 
two phases divided by the BKT critical line. They find that, where the Z symmetry is 
unbroken, the system is in a confining phase, more precisely in a "Coulomb phase", where 
electric charges mutually interact by a logarithmic Coulomb potential. This logarithmic 
behavior has as a consequence a power law dependence of the Polyakov loop correlators. 
Indeed, if the interaction potential is of the form V{r) = alogr, then one can write 
immediately 

G{r) = (n(0)n*(r)) = e"^ = e"^^"^*^ = r"^ = r'^^^) . (14) 
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At the tree level, they find 

G(r)trec oc r~3fef . (15) 
From Eqs. ( fT^ and ( fTSl) the value of r] can be easily found: 

V = ■ (16) 
' 47rT ^ ' 

Above Tc the electric charge is not confined and the system is in a Debye plasma phase. 
The expected large distance behavior for the connected Polyakov loop correlators in this 
case is given by 

Gconn(r) OC 6"*^^^ . (17) 

The Debye mass makes the Coulomb interaction short-ranged and fermions and 
antifermions are approximately free particles. 



3 The lattice QED3 theory 

In this paper we discretize the Euclidean action on a lattice with N'^ x A^,- sites and lattice 
spacing a by staggered fermions X) according to 

N 

^ = 5^ + ^^^^(n)M„,„x.M , (18) 

1=1 n,m 

where Sq is the gauge field action and the fermion matrix is given by 

Mr,,m[U]= ^^{[^-(^)]^™,n+.-[f^i(^)]^™,n~p} ' (19) 

v=l,2,3 

with ri,y{n) = [—1)^p<''^p and U^(n) = e*"''*^'^) is the link variable; the phase a is related 
to the gauge field by a^{n) = agA^{n). 

In a non-compact formulation of QED3, Sg is given by 

SG\a] = f E FAn)F,An) , (20) 

where 

F^,u{n) = {a^{n + fi) - a^{n)} - {a^{n + 0) - a^(ra)} (21) 

and afj_{n) is the phase of the link variable and f3 = l/{g'^a). 

With this action we simulate A^ = 1 fiavors of staggered fermions corresponding to 
Nf = 2 fiavors of parity-invariant four- component fermions. 
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3.1 Order parameters and other lattice operators 

The lattice version of the operator ([H]) is ^2] 



UUx) = lle^'^^'-\ (22) 



T=l 



where 6 = gAo and x hves on the spatial lattice, from which we can determine the 
operator averaged on the lattice configuration, 

n- = 4iEn™(^)- (23) 



X 



As noted in Ref. [12], the original Z symmetry is translated into a symmetry on the 
lattice operator 11^; in the broken phase, the operator Um will fluctuate around the values 
e*~ , where A; = 0, 1, . . . , m — 1. In order to study the breaking of this symmetry it is 
useful therefore to introduce the following operator: 

(\ m 
^X]n„(f)j . (24) 

Another possible choice of lattice order parameter is [l3] 

6 = cos [m X arg (11^)] ; (25) 

both definitions (|2^ and (|25ll have the effect to "rotate" the m different phases to the 
direction corresponding to ^ = 0. 

There are two other quantities of physical relevance to be considered; one is the 
monopole density p, which can be defined in the same way as in the compact theory, 
using the method of Ref. jll] . 

P=2 NIK ' ^''^ 

where Nm (A^m) is the number of monopoles (antimonopoles); the other is chiral conden- 
sate, (xx)- 

The presence of transitions is detected by looking for peaks in the susceptibility of the 
operators fl^ . fl2^ and fl^ . the susceptibility of a generic operator O being defined as 

Xo = {O') - {Of . (27) 

The susceptibilities of monopole density and chiral condensate is also considered in order 
to study the possible relation of these latter operators with the confinement / deconfinement 
transition. 
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3.2 General strategy of the lattice calculation 

The main goal of the numerical computations in this paper is to find evidences of the 
existence of a transition between two phases, one with unbroken and the other with 
broken Z symmetry. The strategy for that is to scan the temperature for a fixed value of 
the bare fermion mass aM and on lattices with given extension to study the behavior of 
the lattice order parameters defined in the previous subsection and of their susceptibility. 
On the lattice, we have 

^'Wa'W/' '''' 
therefore at fixed coupling constant g and on a lattice with fixed N^, the temperature can 
be changed by changing j3 (note that the theory is super-renormalizable, therefore g does 
not depend on a). We expect that changing (3 we meet a clear signature of the transition, 
in correspondence to the critical temperature, which, at infinite fermion mass, is given in 
Eq. (|T2l) and, through (1281) . can be translated on the lattice to the following critical value 
for the (3 parameter: 

/3,(M ^ oo) = . (29) 

OTT 

The next step is to study the behavior of the Polyakov loop correlator. If the theoretical 
expectations introduced in Section 12.21 are true, then for a (3 value below the critical 
one, f3c{M), we should find a power law behavior, whereas above (3c{M) there should be 
exponential fall-off. In the Coulomb-confined phase, i.e. for (3 < (3c{M), where the power 
law behavior is expected, we can compare our findings with the expected tree-level value 
of 7] given in Eq. f|T6l) . which in lattice units reads 

N 

V = ■ (30) 

We recall that 1/m represents the ratio g/g., i.e. the fraction charge in units of the 
fundamental charge. It should be noted here that rj is proportional to 1/m?. 



4 Numerical results 

In this Section we present our numerical results obtained on lattices 12^ x 8, 32^ x 8 
and 64^ x 8. On the smallest lattice considered, we found that tunneling effects among 
the different vacua of the theory in the broken symmetry phase are evident in the 
expectation value of Tim- For this observable this makes possible to better unravel the 
vacuum structure in the deconfined phase. However, for the study of correlation functions, 
where tunneling in the deconfined phase is potentially dangerous, we used lattices with 
= 32 and 64 for which we verified that tunneling is absent. 
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4.1 The algorithm 

We used the Hybrid Monte Carlo Algorithm built by superimposing the Metropolis 
acceptance test to the Hybrid Algorithm in the structure proposed in Ref. [IE]). Here is 
a list of the internal parameters of the Monte Carlo code: 

• the mass term uj for the pseudofermionic field in the Hamiltonian, which drives the 
molecular dynamics; 

• the integration step 5t, fixed to 0.03 in all simulations; 

• the frequency MCR of the algorithm refreshments and Metropolis tests. 

We have verified the important role of uo in simulations: if uo is too small {uj <^ 1), the 
system thermalizes very slowly, but observables do not fluctuate much around the mean 
value during the simulation time; on the contrary, if u is not so small (a; < 1), the system 
thermalizes very quickly with the disadvantage that fluctuations increase considerably. It 
is therefore fundamental to optimize the choice of the parameter u. 

Additionally, we use the following abbreviations: FOM for the frequency of measure- 
ments and NMS for the number of measurements or statistics. The bare fermion mass is 
fixed to aM = 0.05. 

4.2 Order parameter 

In Figure [T] we show scatter plots of the order parameter H^ in the complex plane for a 
f3 value known a posteriori to lie in the confined phase and for several values of m. For 
m = 1, 2 data distribute on a disk, while for m > 3 they are on a ring whose radius 
increases with m; in both cases data are uniformly spread around the origin. Figure [2] 
gives the evolution in the simulation time of Hio; it is interesting to observe that argHio 
spans in a continuous way the interval [0, 2tt] (with account of the 2n periodicity). 

Figure [3] is the same as Figure [H, but for a value of f3 in the deconfined phase. In 
this case the symmetry is broken and the values of the order parameter accumulate 
in correspondence of the m roots of the identity. Accordingly, the evolution in simulation 
time of argHio spans the interval [0, 27r] in jerky way (see Figure HI). 

These scatter plots tell us that real and imaginary parts of (Hm) are zero in both 
phases, if, evidently, there is enough tunneling in the deconfined phase and if the algorithm 
explores the whole configuration space in the confined phase, as it seems to be the case 
on the lattice 12^ x 8 used to obtain the above Figures. So, the more informative quantity 
here seems to be (abs(Hm)). The behavior of this quantity is shown in FigureOfor varying 
(3 and for several values of m and in Figure [6] for varying m and for different values of f3. 
Figure shows that there is a smooth increase with f3, more evident for higher values of 
m. Data do not allow to clearly single out any transition; indeed, also the susceptibility 
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of abs{Ilm) does not show clear peaks for varying (3 in the same region (see Figure [7]). 
Therefore, we can conclude that this operator is not convenient to distinguish the two 
phases which lead to the different behaviors shown in Figures [T] and [3l 

This conclusion suggests to move to the operator 11™, which evidently has no effect 
in the angular distribution in the confined phase, i.e. in the situation of Figure [1], while 
"concentrates" the angular distribution along the first of the m roots of the identity in 
the deconfined phase, i.e. in the situation of Figure El Equivalently, we can say that 
(argll^) = in the deconfined phase. The operator IIJ^ should therefore permit to 
distinguish between broken and unbroken symmetry phases; since Im(n™) is always zero, 
we consider from the beginning Re(n5^). 

Results are shown in Figures [H] and M In this case we have a clear hint on the position 
of the transition point between the two phases. In particular. Figure M indicates that, 
for all values of m > 1, there is a sharp increase of the signal at (3 ~ 0.3, while Figure [9] 
suggests that the increase is more pronounced the higher is m. For m = 1 data show 
no transition at all, as it must be, since m = 1 means integer electric charge. This is 
confirmed in Figure [10] where a peak in the susceptibility of the operator Rell™ appears 
for all values of m > 1 at a /3c value which, on this lattice, has a mild dependence on m, 
namely, it decreases with m and seems to stabilize for high m around f3c = 0.32. 

In Figures [TT] and [T2] we show the same results obtained on a lattice 32^ X 8. In this 
case the values of (3 where the peaks appear are more stable with m and we can estimate 
the critical value as Pc = 0.33. This leads us to argue that the m-dependence of Pc seen 
on the smaller lattice could be a finite volume effect. 

It is important to compare the numerical result for j3c with the theoretical expectation 
discussed in Section 12. 2^ in particular with the value of (3c given in Eq. [22 (31^ = tt^^ ~ 
0.318. . .. This value is impressively close to our numerical result, thus supporting the 
theoretical scenario proposed in Ref. [2] and the conjectured BKT transition (see also 
Ref. [n]). 

In Figures [12] and [Tl]we consider the observable (cos[m x argils]), which should give 
similar information as the previous observable, concerning the location of the transition 
point. We see that data for m > 1 fall on top of each other, while data for m = 1 describe 
a different curve, but also exhibit a sharp increase at the same point as data for m > 1. 
Since we know that for m = 1 the observable cannot be an order parameter, this effect 
should be an artifact of the finite volume. Unfortunately, the susceptibility, shown in 
Figure [T5l does not have a peak structure, but shows rather a jump between two constant 
values, therefore it does not allow for an accurate determination of (3c- 
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4.3 Polyakov loop correlators 

We have studied the wall-wall correlation between two Polyakov 



G(r) = Re(n(0)n*(r)) 



(31) 



for two values of (3, (3 = 0.25 and (3 = 0.40, the first below the critical temperature and 
the other above it. 

It is worth noting here that the dimensionless parameter g'^/M is equal to 80 for 
(3 = 0.25 and to 50 for (3 = 0.40. The other important dimensionless parameter T/M is 
equal to 2.5 in all simulations performed in this work. These values for g'^/M and T/M 
are far from the large M limit, where the theoretical analysis of Refs. [2], [3] was carried on 
(see Subsection 12. 2p . However, the authors of Ref. [3] state that they expect their results 
are stable toward lower electron masses. 

For j3 = 0.40 we have simulated the system on two different lattices: 32^ x 8 with 
uj = 0.2, MCR=10, FOM=10, NMS=100000; 64^ x 8 with u = 0.1, MCR=50, FOM=250, 
NMS=8000. Data for the correlator, shown in Figure [16] (see also Figure[T7]for the scatter 
plot of Um) have been fitted, in a range [rmin,rmax], using the law 



for both lattices and for the different values of m we have obtained x^/d.o.f. < 1. We 
have done "sliding window" fits, that is we have varied the value of Vmin and Vmax until 
obtaining stable values for the fit parameters. 

The most important result is that aAi depends neither on the value of m, and so 
on the fractional charge, nor on the lattice size: aAi = 0.20 ± 0.05. It is instructive to 
determine the "correlation length" related to this mass: ^/a = l/{aAi) = 5.00 ± 1.25; 
this value, much smaller than the lattice size ^/a <^ N„, justifies a posteriori the fit with 
an exponential function. We can conclude that this value has a physical meaning: it is 
the Debye mass characterizing the phase above the critical temperature. 

Also for f3 = 0.25 we have studied the system on two lattices: 32^ x 8 with u = 0.1, 
MCR=50, FOM=1000, NMS=10000, and 64^ x 8 with u = 0.1, MCR=50, FOM=250, 
NMS=10000. 

As discussed in Sec. 12. 2[ we expect that here the best fit to the correlator should be 
given by a power law (see Eq. ( JT4l) ): in any case, to be conservative, we consider both the 
exponential and the power law. Therefore we should interpolate data shown in Figure fTSl 
with the following two functions: 



G(r) = A (e- 



+ e 



-MiN^-r)\ _^ Q . 



(32) 



G(r) 



G(r) 



A {e~^- + e-^(^-'-)) + C , 
: A {r-" + {N, - r)-'') + C . 



(33) 



(34) 



^Note that the value of G(r) = lm(n(0)n*(r)) is always compatible with zero. 
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Here the constant C should zero, since the Z symmetry is unbroken, however, we included 
it to take into account the possibility that, due to not enough large statistics, the whole 
configuration space is not explored by the simulation algorithm, this leading to (Jim) 7^ 
(see Figure fTOil . We note, however, that the large distance correlator does not go to zero 
(see Figure fTSll : this happens because of the small values of M. and rj. In this case, if 
we try to fit data using directly Eqs. flH^ or flMl) . we find that the value of C is not 
equal to the expected value, i.e. abs((nm))^. Therefore, we first determine the connected 
correlator, by subtraction of the numerical value of abs((nm))^, and then fit data using 
Eqs. ( l33l) or ( IMl) without the additive constant C: 



In Figure [20] we show the value of r] obtained on the lattice with N„ = 64 by a fit with 
Eq. flHUl) for different values of m. Simulations on the lattice with = 32 do not permit 
to obtain reliable estimates, since results are not stable. 

It is interesting to note that the behavior of r] is the one predicted by Eq. ( l30l) . except 
for an additive constant, since we find 77 = A/m"^ + B (see Figure I^Tj) : for the lattice with 
N„ = 64, we find A = 1.59(64) and B = 0.270(27), with xVd.o.f. ~ 0.01. Note that the 
tree level theoretical value is 1] = Nr/ (m^47r/9) = 2.5464 . . . /m^. 

If we try to fit data with Eq. (1351) . for both lattices the result is a value of aA4 that 
depends weakly from m and it is of the order 1/N^\ this is a strong evidence that the 
system is in a critical region with a "correlation length" ^/a oc N^j, that is for < T < Tc 
the system is always critical. This is a typical feature of the BKT transition. 

The error bars on aM. and 77 are not those resulting from the fit, which would be 
underestimated owing to the correlation among the values of G{r) at different r's. They 
are estimated, instead, through the behavior of the effective rj and aA4, built from suitable 
ratios of G{r) in such a way that the dependence on the parameters A and B disappears. 
This method allowed also to cross-check the results of the fit. 

4.4 Chiral condensate 

The results on the chiral condensate are presented in Figures [22] and [23 Data show a 
smooth transition for increasing f3, which supports the argument that there is a chiral 
transition coinciding with the confinement-deconfinement one. Results in favor of this 
conclusion can be found also in Refs. [IHl HH] • 

One would expect, however, that the presence of a transition should be accompanied 
by a peak in the susceptibility. This seems to be not the case in our analysis - see 
Figure [23] A possible reason for this unexpected result could be the small volume used. 



G(r) = A(e-^'' + e-^(^'^-'^)) , 



(35) 



G(r) = A {r-" + {N^ - ry) . 



(36) 
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4.5 Monopole density 

The behavior of the monopole density and of its susceptibihty with (3 is shown in Figures!^ 
and [25], respectively. The result is somewhat surprising: the monopole density decreases 
very rapidly in the same region where the Polyakov loop operators show a fast change, 
but the peak in the susceptibility is located at a different value of the coupling constant, 
i.e. at pp ~ 0.18. 

5 Conclusions and outlook 

In this paper we have presented an extended numerical analysis on the confinement- 
deconfinement transition in non-compact QED3 with massive fermions at finite tempera- 
ture. 

We have studied the system for a given value of the fermion mass on lattices with 
different extensions. We have found compelling evidence that there is indeed a transition 
temperature from a high-temperature phase where fractional charges are deconfined or 
Debye-screened and a low-temperature phase with Coulomb-confinement. To detect the 
transition point we have adopted some suitable lattice operators, sensitive to the breaking 
of the underlying symmetry of the system, the Z symmetry. The wall-wall correlation of 
these operators has permitted to characterize the two phases: the confined one exhibits 
power law fall-off with the distance, whereas the deconfined one shows exponential decay. 

There are several indications that the transition and the phase structure are compatible 
with the suggestions of Refs. [21 E] : 

• the critical temperature found by numerical simulations is in remarkable agreement 
with the one estimated for large fermion mass in Ref. [2] and, in its turn, consistent 
with the BKT scenario, = g"^ / (Svr); this agreement may be either an indication of 
smooth mass dependence of the effective action used in Ref. [2] or an accidental fact; 
in both systematic study of the dependence of our results on the fermionic 
mass should be performed; this is, however, beyond the scope of the present work; 

• two phases are clearly seen, one where the correlator of the order parameter exhibits 
a fall-off with the distance with a power law, the other where the fall-off is expo- 
nential; the fact that the power law is valid well inside the confined phase and not 
only on the critical point is an indication in favor of the BKT scenario; 

• the scaling of the index rj for the power law fall-off with the fractional electric charge 
is in agreement with the prediction from Ref. [3j (except for an additive constant), 
which supports the BKT scenario. 

There is an indication that the chiral transition has a relation with the deconfinement 
transition, which should, however, be confirmed by an analysis on larger lattices. 
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If one defines monopoles on the lattice in the non-compact theory as in the compact 
theory, one can see that their density across the deconfinement transition shows a sharp 
change. Surprisingly enough, the susceptibility of the monopole density shows a peak at 
a smaller temperature than the deconfinement one. 

The future development of this work includes a finite-size scaling analysis near the 
transition temperature, in order to achieve an accurate enough determination of the criti- 
cal indices and to conclude that the transition is definitely BKT. Moreover, simulations for 
different values of the fermion mass could allow to explore the phase diagram of the theory 
in the space of the mass, temperature and coupling parameters and to get a unifying view 
of the statistical properties of the system. 
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Figure 1: Scatter plots of 11^ (Imlljn versus ReHj^) for 8 values of m on a 12^ x 8 lattice at 
P — 0.15 (confined phase); the simulation parameters are aM — 0.05, u; — 0.1, St — 0.03, 
MCR=50, FOM=1000, NMS=10000. For m = 1,2 data are homogeneously distributed 
on a disc; for m > 3 data are homogeneously distributed on a ring whose radius increases 
with m. 
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Figure 2: Evolution with Monte Carlo time of modulus (down) and argument (up) of IIio 
on a 12^ X 8 lattice at /3 = 0.05 (confined phase); the simulation parameters are as in 
Figure [H 
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Figure 3: Scatter plots of (Imllm versus Rellm) for 8 values of m on a 12^ X 8 lattice 
at /9 = 0.80 (deconfined phase); the simulation parameters are aM = 0.05, u = 0.1, 
6t = 0.03, MCR=50, FOM=1000, NMS=10000. In this case data are distributed on the 
m roots of unity. 
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Figure 4: Evolution with Monte Carlo time of modulus (down) and argument (up) of IIio 
on a 12^ X 8 lattice at (3 = 0.80 (deconfined phase); the simulation parameters are as in 
Figure [31 



1.0 
0.9 
0.8 
0.7 
"sO.6 

0.3 
0.2 
0.1 
O.Q 



I I I I I I I I I I I I I I I I I I I I I i_i I j_ I I I I I I ^ I I I I I I I 



I ,y J, I ' 




I I I I I I I I I I ■ I I I I I I I I ■ I I I I I I I I I I I I 




.0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 



Figure 5: Modulus of (11^) for 10 values of m versus /3 on a 12^ x 8 lattice; the solid 
line is to guide the eye; the simulation parameters are aM = 0.05, u = 0.1, 6t = 0.03, 



MCR=50, FOM=1000, NMS=10000. 
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Figure 6: Modulus of (11^) for various values of (3 versus m on a 12^ x 8 lattice; the solid 
line is to guide the eye; the simulation parameters are as in Figure O 




Figure 7: Susceptibility of absfllm] for 10 values of m versus /3 on a 12^ x 8 lattice; the 
solid line is to guide the eye; the simulation parameters are as in Figure [5l 
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Figure 8: Real part of (H™) for 10 values of m versus /3 on a 12^ x 8 lattice; the solid 
line is to guide the eye; the simulation parameters are aM = 0.05, u = 0.1, 6t = 0.03, 
MCR=50, FOM=1000, NMS=10000. 
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Figure 9: Real part of (H™) for various values of f3 versus m on a 12^ x 8 lattice; the solid 
line is to guide the eye; the simulation parameters are as in Figure [HI 
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Figure 10: Susceptibility of Re[nj^] for 10 values of m versus /3 on a 12^ X 8 lattice; the 
solid line is to guide the eye; the simulation parameters are as in Figure [HI 
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Figure 11: Real part of (H^) for 8 values of m versus /3 on a 32^ x 8 lattice; the solid 
line is to guide the eye; the simulation parameters are aM = 0.05, u = 0.2, 6t = 0.03, 
MCR=10, FOM=10, NMS=100000. 



22 



0.04 




"•"O.O 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 

Figure 12: Susceptibility of Re[n™] for 8 values of m versus /5 on a 32^ x 8 lattice; the 
solid line is to guide the eye; the simulation parameters are as in Figure [HJ 
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Figure 13: (cos[m x arg 11^]) for 10 values of m versus /3 on a 12^ x 8 lattice; the simulation 
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Figure 14: (cos[m x argils]) for various values of P versus m on a 12^ X 8 lattice; the 
solid line is to guide the eye; the simulation parameters are as in Figure [ 
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Figure 17: Scatter plots of 11^ for 10 values of m on on a 64^ x 8 lattice at /5 = 0.40 
(deconfined phase); the simulation parameters are aM = 0.05, u = 0.1, 6t = 0.03, 
MCR=50, FOM=250, NMS=8000. 
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Figure 18: G{r) versus r on a 64^ x 8 lattice at /? = 0.25 (confined phase); the simu- 
lation parameters are m = 10, aM = 0.05, = 0.1, 6t = 0.03, MCR=50, FOM=250, 
NMS=10000. 
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Figure 19: Scatter plots of 11^ for 10 values of m on on a 642 X g lattice at /3 = 0.25 
(confined phase); the simulation parameters are aM = 0.05, u = 0.1, 6t = 0.03, MCR=50, 
FOM=250, NMS=10000. 
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Figure 20: rj versus m on a 64^ x 8 lattice at (3 = 0.25 (confined phase); the simulation 
parameters are as in Figure | 
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Figure 21: 1] versus on a 64^ x 8 lattice at /5 = 0.25 (confined phase); the simulation 
parameters are as in Figure [181 
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Figure 22: Chiral condensate versus /3 on a 12^ x 8 lattice; the simulation parameters are 
aM = 0.05, uj = 0.1, 6t = 0.03, MCR=50, FOM=1000, NMS=10000. 
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Figure 23: Chiral susceptibility versus /9 on a 12^ x 8 lattice; the simulation parameters 
are as in Figure! 
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Figure 24: Monopole density versus P on a 12^ x 8 lattice; the simulation parameters are 
aM = 0.05, u = 0.1, 6t = 0.03, MCR=50, FOM=1000, NMS=10000. 
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Figure 25: Monopole density susceptibility versus /3 on a 12^ x 8 lattice; the simulation 
parameters are as in Figure | 



30 



